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Abstract. We calculate a grid of models with and without the effects of axial rotation for massive stars in the range 
of 9 to 60 Mq and metallicity Z = 0.004 appropriate for the SMC. Remarkably, the ratios fi/fi C rit of the angular 
velocity to the break-up angular velocity grow strongly during the evolution of high mass stars, contrary to the 
situation at Z = 0.020. The reason is that at low Z, mass loss is smaller and the removal of angular momentum 
during evolution much weaker, also there is an efficient outward transport of angular momentum by meridional 
circulation. Thus, a much larger fraction of the stars at lower Z reach break-up velocities and rotation may thus 
be a dominant effect at low Z. The models with rotation well account for the long standing problem of the large 
numbers of red supergiants observed in low Z galaxies, while current models with mass loss were predicting no 
red supergiants. We discuss in detail the physical effects of rotation which favour a redwards evolution in the HR 
diagram. The models also predict large N enrichments during the evolution of high mass stars. The predicted 
relative N-enrichments are larger at Z lower than solar and this is in very good agreement with the observations 
for A-type supergiants in the SMC. 
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1. Introduction 

The aim of this work is to study the effects of rotation in 
massive stars at low metallicity Z as in the SMC. There 
are 3 main problems emerging in this context. 1.— Firstly, 
we want to examine whether there are differences between 
the effects of rotation at low Z and at solar metallicity. 
Surprisingly, we find some striking differences, which sug- 
gest that rotation may be a much more important ef- 
fect in the evolution of stars at Z lower than solar. 2.- 
Then, we want to address the problem of the number ra- 
tio B/R of blue to red supergiants, which is one of the 
most severe problems in stellar evolution. As stated by 



Kippenhahn and Weigert (1990) about stars on the blue 
loops in the He-burning phase, the present phase " ... is 
a sort of magnifying glass, revealing relentlessly the faults 
of calculations of earlier phases." As long as the problem 
is not solved, we cannot correctly predict the populations 
of massive stars in galaxies at various metallicities Z, nei- 
ther the spectral evolution of galaxies, nor the chemical 
yields. 3.- A third important question concerns the chem- 
ical abundances of massive stars at low metallicities. In 
particular, we recall that very large nitrogen enrichments 
have been found for the A-type supergiants in the SMC 



(Venn 1998 1999 ). There, the relative N/H excesses can 
reach a factor of 10, much larger than in the Galaxy. This 
is a strong constraint, which needs to be examined care- 
fully. 

Let us comment a bit more on the embarassing prob- 
lem of the number rati o B/ R, which has been reviewed 
by L an ger & Maeder ( 1995| ) and by Maeder & Meynet 
( 2000a ). To make a long story short, the observations of 
nearby galaxies show that there are much more red super- 
giants at metallicities lower than solar, so that B/R de- 
creases steeply with decreasing metallicity (Humphreys & 
McElroy 1984 ) . Large differences in the B /R ratio also ex- 
ist between clusters in the Galaxy and in the SMC, where 
the red supergiants were found more numerous by an order 
of a magnitude (Meylan & Maeder 1982). 

There are no sets of models which correctly predict the 
observed trend of decreasing B/R with decreasing Z, as 
emphasized by Langer & Maeder (1995). All kind of mod- 
els were examined: models with Schwarzschild's criterion 



(cf. Stothers & Chin |1992[ ), models with Sc hwarz schild's 
criterion and overshooting (cf. Schaller et al. 1992 ), mod- 
els with Ledoux crit erion (cf. Stothers & Chin 1992 ; 



Brocato & Caste llani 1993 ), models with semiconvection 
(cf. Arnett 1991), models with semiconvective diffusion 
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(cf. Langer & Maeder 1995). The comparisons generally 
show that the models with Schwarzschild's criterion (with 
or without overshooting) may reproduce the observed 
B/R at solar metallicity, while they fail at lower Z. At 
the other extreme, the models with the Ledoux criterion 
and those with semiconvection reproduce well B / R at the 
metallicities of the SMC, but they fail at higher Z. This 
shows that, at least in part, the problem is related to the 
size of the core and to the mixing efficiency outside the 
core. 

Of course, other effects such as mass loss, convection, 
opacities and metallicities play a role. In particular, as 
evide nced by the Ge neva g rids of models (Meynet et al. 



used here, with a few improvements, which are mentioned 
below. 



2.1. Meridional circulation, shears and horizontal 
turbulence 

Meridional circulation plays a major role in the redistribu- 
tion of the angular momentum in stars. Let us write the 
equations for the radial term of the vertical component 
U(r) of the meridional circulation (cf. Maeder & Zahn 
|1998 ), which will be essential for the following discussions, 



1994| ; cf. also Maeder |1981| ) a growth of mass loss favours U(r) = 



the formation of more red supergiants in the He-burning 
phase. However, the B/R problem in the SMC cannot be 
solved by mass loss, because the mass loss rates at lower 
Z, as in the SMC, are smaller than in the Galaxy and this 
produces fewer red supergiants. 

In some recent works, it has been shown that the ac- 
count of the axial rotation of stars changes all the model 
outputs; tracks in the HR diagram, lifetimes, surface abun- 
dances, etc... (Meynet & Maeder 2000| ; Heger & Langer 
200C ) . Noticeably rotation, both by its effects on internal 
mixing and on mass loss was found to favour the redward 
motions. This is a positive indication and we now further 
explore it. A grid of stellar models with rotation in the 
range of 9 to 60 M Q and metallicity Z — 0.004 was con- 
structed. It will provide a basis for comparison with the 
SMC observations. 

In Sect. 2, we briefly discuss the improvements brought 
in the model physics. In Sect. 3, we discuss the evolution of 
the internal rotation law f2(r) and in Sect. 4 of the surface 
rotation velocities v. The models with zero rotation are 
briefly discussed in Sect. 5. In Sect. 6, the results for the 
HR diagram and the lifetimes are shown. Sect. 7 is devoted 
to the study of the B/R ratio from the numerical models; 
we also account for the results in terms of the physical 
properties of stellar models. In Sect. 8, we analyze the 
chemical abundances at the stellar surface and compare 
the results to observations of A-type supergiants in the 
SMC. Sect. 9 gives the conclusions. 

2. Physics of the models 

We apply the treatment of the hydrostatic effects appro- 
priate to differential rotation, as described by Meynet & 
Maeder (1997). We recall that the classical treatment fre- 
quently applied is as a matter of fact not correct in the 
presence of shellular differential rotation. We also account 
for the rotational distortion and for the von Zeipel theo- 
rem. The T e g given here corresponds to the value as ob- 
served for an average orientation angle between the axis 
of rotation and the direction of the observer. We include 



P 



the effects of shear diffusion (Zahn [199% Maeder [1997D 
and the effects of the meridional circulation, as studied 



by Maeder & Zahn ( [19981) . The same p hysics of the mod- 
els as described in Meynet & Maeder ( |2000 , paper V) is 



WCpT [V ad - V 
L 

m7 



'" J+ T¥ 



(1) 



with the current notations as given in the quoted paper. 
The main term in the braces in the second member is Eq. 
If we ignore secondary terms, it behaves essentially like, 
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The overlined expressions like p mean the average over the 
considered equipotential. The term with the minus sign 
in the square bracket is the Gratton-Opik term, which 
becomes important in the outer layers due to the decrease 
of the local density; it can produce negative values of U(r). 
A negative U (r) means a circulation going down along the 
polar axis and up in the equatorial plane, thus making an 
outwards transport of angular momentum. The positive 
or negative values of U(r) play a major role in massive 
star evolution. 

In paper V we did not account for the effects of hori- 
zontal turbulence on the shears. This was done by Talon 



& Zahn (1997). They found that the diffusion coefficient 
for the shears is modified by the horizontal turbulence. 
The change can be an increase or a decrease of the diffu- 
sion coefficient depending on the various parameters, as 
discussed below. Thus, we also have to include the devel- 



opments by Talon & Zahn (1997), we have 
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with D h ~ \rU{r)\ 



(4) 



wher e is the coefficient of horizontal diffusion (cf . Zahn 
1992 ). We ignore here t he thermal coupling effects dis- 
cussed by Maeder ( 1997 ) because they were found to be 
relatively small and they increase the numerical complex- 
ity. Interestingly, we see that in regions where V M ~ 0, 
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Eq. (3) leads us to replace K by (K + D\/) in the usual ex- 
pression (cf. Talon & Zahn 1997), i.e. it reinforces slightly 
the diffusion in regions which are close to chemical homo- 
geneity. On the contrary, in regions where V M dominates 
with respect to (V a< a — V ra d), the transport is proportional 
to Dh rather than to K, which is quite logical since the 
diffusion is then determined by rather than by ther- 
mal effects. The above result shows the importance of the 
treatment for the meridional circulation, since in turn it 
determines the size of Dh and to some extent the diffusion 
by shears. The equation of transport of angular momen- 
tum and of chemical elements are followed explicitly. 

2.2. Mass loss rates and their dependence on rotation 

For t he ma ss loss rates, the recent data by Kudritzki and 
Puis ( EOOOl) for OB stars and B-type supergiants are ap- 
plied, they are completed f or A- F supergiants by the ex- 
pressions by de Jager et al. ( 1988 ). For the red supergiants, 
there are a numb er of recent parametrizations proposed 
(cf. Willson |2000|) , from low to high rates. We choose fi- 
nally the medium -high rates as they were proposed by 
Vanbevere n et al . ( 1998 ). For WR stars, the rates by Nugis 
& Lamers ( 2000 ) were applied. These rates account for the 
clumping of matter in the winds from WR stars and they 
are lower by about a factor of 2 tha n those used in the 
Geneva grids (cf. Schaller et al. 1992 ). 

We have to account for rotation effects on the mass loss 



rates. In our previous work (Meynet fc Mae der 200C ), we 
used the expression derived by Langer ( 1998 ), based on the 
numerical models by Friend & Abbott (198E). However, 
these models did not include the gravity darkening pre- 
dicted by the von Zeipel formula. A new expression has 
been derived by the application of the wind theory over 
the surface of a rotating star, taking also into account var- 
ious improvements in the model of a rotating star (Maeder 
& Meynet 2000b ). The ratio of the mass loss rate of a star 
rotating with an angular velocity O to that of a non rotat- 
ing star at the same location in the HR diagram behaves 
as 
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The term T is the usual Eddington factor for electron scat- 
tering opacities and p m is the average density internal to 

o 2 

the surface equipotential. The term 2 kG p ls a function of 
the ratio of the rotational velocity to the break-up velocity 



fi 2 



2nGp n 



(6) 



0.8. Here, 



with deviations of less than ~ 3% up to — 

(O /— y ji j \ 2 
3 7j-^ ) and R p b is the polar radius at break-up 

(i? p b can be taken equal to the radius of the non-rotating 
star). The coefficient a in Eq. (5) is a force multiplier, 



which depends on T e g. We use the empirical force multi- 
pliers derived by Lamers et al. fll995| ): a = 0.52, 0.24, 0.17 
and 0.15, for logoff > 4.35, 4.30, 4.00 and 3.90 respec- 
tively. 

Of course, we have to account for the fact that the em- 
pirical values for the mass loss rates used for non-rotating 
stars are based on stars covering the whole range of rota- 
tional velocities. Thus, we must apply a reduction factor 
to the empirical rates to make them correspond to the non 
rotating case. After convolution of the effects described by 
Eq. (5) over the observed distribution of rotational veloc- 
ities, taking also into account that the axes of orientation 
are randomly distributed, we estimate that the average 
correcting factor is equal to about 0.8, which has to be ap- 
plied to the mass loss rates for the Main-Sequence (MS) 
OB stars. For post-MS stars, the mass loss rates are too 
uncertain anyway, and the multiplying correction factor is 
almost unity, so that we do not apply any correction. 

2.3. Initial composition, opacities, nuclear reactions 

For a given metallicity Z (in mass fraction), we esti- 
mate the initial helium mass fraction Y from the relation 
Y = Y p ■ AY/AZ Z, where Y p is the primordial helium 
abundance and AY/AZ the slope of the helium-to-metal 
enrichment law. Values of Y p between 0.23 and .24 are 
given in th e literature (see e.g. Pagel et al. 1992 ; Izotov 
et al. |1997| ; Peimbert et al. |2000|) . Observations of HII re- 



gions in blue compact dwarf galaxies indicate values for 
AY/AZ equal t o 1.7 ±0.9 (Izotov et al. |1997| ), 2.3±1.0 
(Thuan & Izotov 1998). Such values are in agreement with 



the AY/AZ value of 2.3±1.5 derived by Fernandes et 
al. (1998) from the study of nearby visual binary stars. 



Here we set Y p = 0.23 and AY/AZ = 2.25 as in the 



recent grids of stellar models by Girardi et al. (2000). 
For the metallicity Z — 0.004 considered in this work, 
we thus obtain X = 0.757 and Y = 0.239. For the mix- 
ture of the heavy elements, we adopted the same mixture 
as the one used to compute the opacity tables. Opacities 



are from Iglesias & Rogers (1996) complemented at low 
temperatures with the molecular opacities of Alexander 



( bttp: / /web. physics. twsu.edu/alex/wwwdra. htm ) . 

The rates for the charged particle reactions are taken 
from the new NACRE compilation (Angulo et al., 1999). 
For the reaction 12 C(a,7) 16 and for the range of tem- 
peratures corresponding to the He-burning phase, the 
NACRE rate is a f actor two higher than the rate of 
Caughlan & Fowler ( 1988j ) and amounts to about 80% 
of the rate a dopted in paper V and taken from Caughlan 
et al. ( 1985 ). Finally, let us recall that as in paper V, we 
adopted the Schwarzschild criterion for convection with- 
out overshooting. 

Thus, with respect to paper V, in addition to the im- 
provements of the physics describing the effects of rotation 
on the mixing of the chemical elements and on the mass 
loss rates, we have updated the mass loss rates (see pre- 
vious section), the initial composition, the opacities and 
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Fig. 1. Evolution of the angular velocity Q, as a function 
of the distance to the center in a 20 M© star with u; n i 
= 300 km s _1 and Z = 0.004. X c is the hydrogen mass 
fraction at the center. The dotted line shows the profile 
when the He-core contracts at the end of the H-burning 
phase. 

the nuclear reaction rates. In order to evaluate the impor- 
tance of the changes which are not linked to the physics of 
rotation, we have computed a non-rotating 20 M Q model 
at solar metallicity with the physics included in this pa- 
per (see Fig. || below). We obtain that the present 20 M Q 
model at solar metallicity follows a very similar track in 
the HR diagram as the one presented in paper V. With 
respect to the model of paper V, the MS lifetime is in- 
creased by about 7%, while the helium-burning phase is 
decreased by the same amount. This results in part from 
the slightly increased initial hydrogen abundance in the 
present model (X=0.705 while in paper V, X = 0.680). 
The mass loss during the MS phase has been significantly 
reduced in the present paper. Indeed, the present 20 M Q 
model at solar metallicity has lost during the MS phase a 
little more than one third of a solar mass, while the mass of 
the corresponding model presented in paper V decreased 
by about one solar mass during the same period. 

3. The evolution of the internal rotation and 
meridional circulation 

The initial convergence of the internal Sl-profile was dis- 
cussed in paper V and here the results are the same. The 
main point was that the Sl-profile converges very quickly 
in about 1 to 2% of the MS lifetime towards an equilib- 
rium profile. The same fast convergence occurs for the in- 
ternal distribution of U(r), the radial term of the vertical 
component of the meridional circulation. Then, during the 



Fig. 2. Evolution of U(r) the radial term of the vertical 
component of the velocity of meridional circulation for the 
same model as in Fig. [l| X c is the hydrogen mass fraction 
at the center. 



MS phase Sl(r) decreases slowly keeping a small degree of 
differential rotation, which plays an essential role in the 
shear effects and the related transport mechanisms. We 
recall here the interesting result shown in paper V that, 
in contradiction to the classical Eddington-Sweet theory, 
the velocity U(r) does not depend very much on the rota- 
tional velocities. 

Fig. [I] shows the evolution of the fi-profile during the 
evolution on the MS of a 20 M Q star with Z = 0.004 
and an initial rotation velocity v- m i = 300 km s . Let us 
point out some differences with respect to the models at 
Z = 0.020. For the same initial rotational velocities 
at a given central H-content X c , the values of f2(r) are 
higher at Z = 0.004. One first reason is rather trivial, i.e. 
because the radius is smaller at lower Z, the same v- ln i 
corresponds to a larger angular velocity SI. Indeed, the 
radius of a 20 M Q star at Z = 0.004 is smaller, being 
equal to 0.8 of the radius of a same mass star at Z = 
0.020. This ratio of 0.8 keeps closely the same during the 
whole MS phase. Another reason for the higher 51 (r) lies 
in the smaller mass loss rates, which lead to smaller losses 
of angular momentum, thus favouring a higher internal 
rotation at the end of the MS phase. 

The internal gradients of fl are higher at lower Z due to 
the higher compactness of the star. This means that the 
outwards transport of chemical elements by shears will 
be favoured at lower Z. However, these larger shears do 
not destroy the f2-gradient, since quite generally shears 
are much less important than meridional circulation for 
transporting the angular momentum and shaping S!(r). 
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Fig. 3. Evolution of the surface equatorial velocity as a 
function of time for stars of different initial masses with 
v ini = 300 km s" 1 and Z = 0.004. 
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Fig. 4. Evolution of the ratio r2/f2 cr ;t of the angular veloc- 
ity to the break-up angular velocity at the stellar surface 
for stars of different masses at Z — 0.004. 



Fig. g shows the evolution of U (r) in the same model. 
U (r) is initially positive in the interior, but progressively 
the fraction of the star where U(r) is negative is grow- 
ing. This is due to the Gratton-Opik term in Eq. (2), 
which favours a negative U(r) in the outer layers, when 
the density decreases. This negative velocity causes an 
outward transport of the angular momentum, as well as 
the shears. Due to the higher density in the envelope at 
Z = 0.004, the Gratton-Opik term is less important and 
the values of \U(r)\ are smaller than at a solar composi- 
tion. However, this is partly compensated by the smaller 
stellar radius which reduces the characteristic time for 
transport. In the model of 20 M Q at Z = 0.004, we 
have the timescales r C j rc ~ R/U = 1.5 • 10 7 yr. and 
Tshear ~ -R 2 /-Dshcar = 8.0 ■ 10 7 yr. This confirms that the 
circulation U(r) is, as for models at Z = 0.02, more im- 
portant than the shears for transporting the angular mo- 
mentum. In addition, we see that r c i rc is of the same order 
as the MS lifetime of 9.7T0 6 yr., which allows a significant 
transport to occur. 

4. The evolution of the surface rotation velocities 

The evolution of the rotational velocities v at the stellar 
surface is shown in Fig. ^. We notice that for all masses 
considered here, v remains almost constant until the end of 
the MS phase. Starting from an initial velocity Vi n i — 300 
km s~ x on the ZAMS, the initial convergence brings v 
to about 250 (± 10 km s _1 ) and these values decrease 
to about 200 km s _1 for stars in the range of 9 to 25 
Mq at the end of the MS phase. This behaviour is very 
different from that found for models with Z = 0.02 (paper 



V), where the values of v for the large masses decreased 
drastically (a value as low as 30 km s _1 was found for 
the 60 M© star at the end of the MS). The reason for the 
decrease of v at Z = 0.02 was the large mass loss rate M, 
while here at Z = 0.004 the M-rates are much smaller 
and the stars retain most of their angular momentum. 
This behaviour is reminiscent of that found in models with 



solid body rotation by Sackmann & Anand (1970) and by 
Langer (|1997j , |1998| ). 

Fig. |^ shows the evolution of the fraction of the critical 
angular velocities fp—', these results are striking and po- 
tentially very important. For the model of 9 M©, this frac- 
tion remains almost constant, but for the higher masses, 
the ratio grows a lot and may even reach 1.0 during 
the overall contraction phase at the end of the MS phase. 
This is the opposite of the results of the models with Z 
= 0.02, where the fraction j^— strongly decreased for the 
largest masses due to their high mass loss, which removed 
a huge amount of angular momentum. In the present mod- 
els at Z = 0.004, mass loss is insufficient to remove enough 
angular momentum. 

We recall that, under the hypothesis of local conserva- 
tion of angular momentum, the ratio jp-r would decrease 
for a growing radius. Thus, the growth of jp— for the 
largest masses in Fig. ^results essentially from the larger 
outwards transport of angular momentum by circulation. 
Indeed, we notice that U(r) in a 40 M Q with Z = 0.004 
near the end of the MS reaches —6 • 10~ 2 cm s _1 , com- 
pared to —0.6 • 10~ 2 cm s _1 in the 20 M© illustrated in 
Fig. H The large negative velocities in the 40 M Q model 
at Z=0.004 transport more angular momentum outwards 
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Fig. 5. Evolution of the surface equatorial velocity as a 
function of time for 20 M Q stars at Z — 0.004 with differ- 
ent initial velocities. 



and thus explain the rise of -^~ t at the surface for large 
masses in Fig. |]. 

Now, why is U (r) more negative in larger masses ? 
From Eq. (1) we see that for given values of (V a d — V ra d) 
and V M , U(r) behaves essentially like 



U(r) 
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(7) 



For increasing stellar masses, we have a smaller value of 
the gravity g, a smaller value of the ratio of P gas to 
the total pressure P, a larger L/M and a smaller density. 
All these differences contribute to make the values of U(r) 
more negative. Together with the smaller mass loss, these 
are the physical reasons for the growth of the surface ratios 
with the stellar masses. These results suggest that, if 
the initial distributions of the rotational velocities do not 
have a higher fraction of slow rotators at lower metallici- 
ties (a point which is still unknown), then we may expect 
that at lower Z, as in the SMC, a much larger fraction 
of O- and B-type stars are close to break-up velocities. 
This suggests that rotation is not necessarily a secondary 
parameter for the evolution of massive stars at low Z. 

Fig. H shows the evolution of the surface rotational 
velocities v as a function of time for different initial ve- 
locities in the case of the model of 20 M at Z = 0.004. 
We notice the rather parallel behaviour of the evolution 
of v. This is rather different from the models at Z — 0.02, 
which showed a certain degree of convergence of v at the 
end of the MS, due to the large mass loss rates. 

With some necessary reservations, let us now make 
a few speculations. For zero metallicity models, the MS 



phase is shifted by 0.27 dex in T e g with respect to the 
models at solar composition, which means radii smaller 
by a factor of 3.4 at a given luminosity for Z = 0. For the 
models at Z = 0.004, the blue shift amounts only to 0.04 
dex, implying a radius smaller by a factor of 1.2. Thus, we 
see that the degree of compactness at very low Z is much 
stronger. Thus, the models at Z = will likely amplify 
the behaviour shown by the models at Z = 0.004, due 
to their higher compactness and their smaller mass loss 
rates. Thus, we may wonder whether a large fraction, if 
not most, massive stars at very low Z do not reach break- 
up velocities during their MS phases. The observati ons o f 
a much higher fraction of Be-stars (cf. Maeder et al. 1999 ) 
in the LMC and SMC compared to the Milky Way tend 
to support this view. 

Finally, we notice in Fig. 4 the occurence of the so- 
called spin-up effect in the 9 M Q on the blue loop: as the 
star is contracting, its rotation strongly accelerates. This 
effect has been described by Heger & Langer ( 1998] ). The 
envelope accelerates its rotation, due to both its contrac- 
tion and the fact that the convective zone is receding, as 
shown by Heger & Langer. 

5. Models with zero rotation 



The evolutionary tracks and the lifetimes of the non- 
rotating tracks at Z — 0.004 are presented in Figs. ^ and]?] 
and in Table |l|. With respect to the previous Geneva grids 
of stellar mo dels a t Z = 0.004 (Charbonnel et al. |1993| 
Meynet et al. [1994 ), those presented here show the follow- 
ing main differences: the MS width is reduced, typically 
at log T c ff = 4.4, the width in AlogL/L© is decreased by 
0.34 dex; for initial mass star models below 40 M , the 
MS lifetimes are shorter by about 11% for the 9 M© model 
and by 1.5% for the 40 M© model; the ratios of the He- 
to the H-burning lifetimes are greater by about a factor 
two for the 9 M Q model and by 7% for the 40 M© model. 
All these effects are quite consistent with the fact that, in 
the present work, we did not include any overshooting (see 
e.g. Bert elli et al. |l985[ Maeder fc Me ynet |l989t Chin & 
Stothers 199C ; Langer & Maeder 1995 for detailed discus- 



sions of the effects of overshooting on the stellar models). 
For initial masses superior to about 60 M©, the effects 
of the stellar winds become very important and dominate 
the effects due to a change of the criterion for convec- 
tion. Qualitatively similar differences are obtained when 
the present results are compared with the stellar models 
at Z = .004 computed with over shoot ing by Fagotto et 
al. ( p^94| ) and Claret & Gimenez ( |l998| ). 

The models with initial masses between 10 and 12.5 
M© are peculiar, they show a behaviour in the HR dia- 
gram intermediate between the cases of stars presenting a 
well developed blue loop (as the 9 and the 10 M© models 
in Fig. Q) and the case of more massive stars which do not 
produce any blue loop, but begin to burn their helium in 
their core at a high effective temperature while they cross 
the HR diagram for the first time (as the 12.5 and 13 M Q 
in Fig. Q). The models in the transition mass range (i.e. 
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the 10.5, 11 and 12 M Q models in Fig. present a "par- 
tial blue loop" which starts at a high effective tempera- 
ture and then the models spend nearly their whole helium 
burning phase in the blue (see Sect. 7). Such a behaviour 
has also been found by Charbonnel et al. (|1993| ) for their 
15 M model and by Claret & Gimenez (1998) for stars 
with initial masses between about 16 and 30 Mq. In this 
last case, the exact mass range depends on the adopted 
value for the initial abundance of helium. These results 
seem to indicate that an extension of the H-burning core, 
e.g. by overshooting, shifts the transition range to higher 
initial masses, which is quite consistent. The models by 
Fagotto et al. ( 1994 ) do not show such a behaviour be- 
tween 15 and 30 Mq. This may be due to the fact that 
these models were computed with a much greater over- 
shooting parameter than the two grids mentioned above. 

6. HR diagram, mass-luminosity relation and 
lifetimes 

6.1. A brief review of the effects of rotation on the 
stellar models 

Many effects induced by rotation at solar metallicity are 
also found here at lower metallicity. T hese effects have 
been des cribed in detail by Heger et al. ( [2000 ), Maeder & 
Meynet ( ^000a| ) and Meynet & Maeder ( p000| ) . Thus we 
will be very brief here. Rotation modifies the evolution- 
ary tracks in the HR diagram through the following main 
physical effects: rotation lowers the effective gravity g e g; 
it enlarges the convective cores and smoothes the chem- 
ical gradients in the radiative zones (the main effect); it 
enhances the mass loss rates; it produces atmospheric dis- 
tortions. Since the star has lost its spherical symmetry, 
the tracks may also change in appearance depending on 
the angle of view. Here we suppose an average angle of 
view, as was done in paper V. 

As a consequence of these effects, depending on the 
degree of mixing (see below), the evolutionary tracks on 
the MS are extended towards lower effective temperatures 
(as would do a moderate overshooting) and/or are made 
over luminous. The evolution towards the red supergiant 
stage is favoured (see Sect. 7), as well as the evolution 
into the WR phase. Moreover, as will be discussed below, 
the MS lifetimes are increased and the surface abundances 
are modified. 

Since there are different rotational velocities, a star 
of given initial mass and metallicity can follow different 
tracks corresponding to various initial rotational veloci- 
ties. To minimize the number of tracks to be computed, 
we have to choose a value of the initial velocity which 
produces an average velocity during the MS not too far 
from the observed value. The problem here is that one 
does not know what this average rotational velocity is for 
OB stars at the metallicity of the Small Magellanic Cloud 
(SMC). In the absence of such an information, we adopt 
here the same initial rotational velocity as in our grid at 
Z = 0.020, namely the value Uj n i = 300 km s _1 . Defining a 



mean equatorial velocity v during the MS as in Meynet & 
Maeder (2000), such a value for v- m i corresponds to values 
of v between 220 and 260 km s _1 . These values are close 
to the average rotational velocities observed for OBV type 
stars at solar metallicity, which are between 200-250 km 



6.2. The HR diagram 

Figure ^| shows the evolutionary tracks of non-rotating 
and rotating stellar models for initial masses between 9 
and 60 Mq. One sees that the MS width is increased, 
as would result from a moderate overshoot. Let us recall 
here that two counteracting effects of rotation affect the 
extension of the MS band (see paper V). On one hand, 
rotational mixing brings fresh H-fuel into the convective 
core, slowing down its decrease in mass during the MS. 
This effect produces a more massive He-core at the end 
of the H-burning phase and this favours the extension 
of the tracks towards lower effective temperatures. On the 
other hand, rotational mixing transports helium and other 
H-burning products (essentially nitrogen) into the radia- 
tive envelope. The He-enrichment lowers the opacity. This 
contributes to the enhancement of the stellar luminosity 
and favours a blueward track. Clearly here, the first effect 
dominates over the second one. 

In this context, it is interesting to recall that, due to 
the account of horizontal turbulence in the present mod- 
els, the mixing of the chemical elements by the shear has 
been effectively reduced in the regions of steep /j,-gradient 
with respect to paper V (see Sect. 2). This favours, for 
a given initial rotational velocity, an extension of the MS 
towards lower effective temperatures. Indeed when rota- 
tional mixing is decreased, either as a consequence of the 
reduction of the initial velocity or by a reduction of the 
shear diffusion coefficient as in the present models, the 
time required for helium mixing in the whole radiative 
envelope is considerably increased, while the time for hy- 
drogen to migrate into the convective core, although it is 
also increased, remains nevertheless relatively small since 
hydrogen just needs to diffuse through a small amount of 
mass to reach the convective core (see Meynet & Maeder 
120001 ). Thus, some increase in the size of the core results, 
while the effect on the helium abundance in the envelope 
is not significant. The same kind of effect can be seen 
in the models by Talon et al. ( |1997| ). Of course the fx- 
gradients are strong near the core and can slow down the 
diffusion process mentioned above, but on the other hand, 
the efficiency of the diffusion of hydrogen will also increase 
with the increasing H-abundance gradient at the border of 
the core. These are the various reasons why the numerical 
models show that, for a moderate rotational mixing, the 
effect of rotation on the convective core mass overcomes 
the effect of helium diffusion in the envelope. 

Let us recall that some core overshooting was needed 
in stellar models in order to reproduce the observed MS 
width (see Maeder & Mermilliod [I98lfc Maeder & Meynet 
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Fig. 6. Evolutionary tracks for non-rotating (dotted lines) and rotating (continuous lines) models for a metallicity Z 
= 0.004. The rotating models have an initial velocity Vi n i of 300 km s -1 . For the rotating 60 M Q model, only part of 
the evolution is plotted. 



1989| ). Typically, the width of the MS band in A logi/L Q 
at log T e ff = 4.4 obtained from the present non-rotating 
models is 0.85, while a value of about 1.20 is required to re- 
produce the observation. Such a width can be reproduced 
with a moderate amount of overshooting (see Schaller et 
al. ~ 



1992) 



Rotation as we just saw above produces also a wider 
MS. Indeed, the width of the MS band in AlogL/L Q at 
log T c ff = 4.4 for the rotating models presented in Fig. || 
is slightly superior to 1. This comparison shows that if 
the ^ini = 300 km s" 1 models are well representative of 
the average case, then rotation alone might account for 
about half of the MS widening required by the observa- 
tions. This is only a rough estimate. The MS extension 
is different for different initial rotational velocities (see 
Fig. 8), it will also be slightly different depending on the 
angle of view of the stars. Therefore the evaluation of the 
contributions of rotation and overshooting to the widen- 
ing of the MS band requires the computations of numerous 
models for various initial masses and rotational velocities 
so that detailed population synthesis models can be per- 



formed. However, we can say with some confidence that 
rotation decreases by about a factor of 2 the amount of 
overshooting needed to reproduce the observed MS width, 
as was already proposed by Talon et al. (1997). 

One striking difference between non-rotating and ro- 
tating models after the MS concerns the fraction of the he- 
lium burning phase spent as a red supergiant. This point 
will be discussed in detail in the next section. Let us also 
mention here that rotation shortens the blue loops of the 
9 M Q model. This is a consequence of the more massive 
helium cores existing at the end of the H-burning phase in 
rotating models, as well as of the addition of helium near 
the H-burning shell (cf. Sect. 7.2). As in the case of the 9 
M Q model and for the same reason, rotation reduces the 
extension of the "partial" blue loop associated with the 
12 M Q model. 

For the rotational velocities considered here, no model 
enters the WR phase during the core H-burning phase. 
But after the MS phase, the rotating 40 M Q model enters 
the WR phase while its non-rotating counterpart does 
not. This illustrates the fact that rotation decreases the 
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Fig. 7. Evolutionary tracks for non-rotating stellar mod- 
els at Z = 0.004 in the transition region between stars with 
and without blue loops. The initial masses are indicated 
in solar masses. 



Fig. 8. Evolutionary tracks for rotating 20 M Q models 
with different initial velocities and various initial metal- 
licities. The initial velocities v; n i are indicated. See Table 
1 for more details on the models at Z = 0.004. 



minimum initial mass for a single star to become a WR 
star (M aeder |1987| Fli egner & Langer |1995| ; Maeder & 
Meynet |2000a|; Meynet |2000[). We shall not develop this 



point further here since it will be the subject of a forth- 
coming paper. 

6.3. Differences between rotating tracks at Z = 0.004 
and 0.020 

Since the physics has been improved with respect to pa- 
per V (see Sect. 2), we cannot directly compare the solar 
tracks of paper V with the present ones. Therefore we 
computed one 20 M model at solar metallicity with and 
without rotation with exactly the same physics as in the 
present paper. These models are plotted together with the 
20 M models at Z = 0.004 in Fig. |. 

From this figure, one sees that rotation at low metal- 
licity has similar effects than at solar metallicity. For in- 
stance, the increase due to rotation of the He-core masses 
at the end of the H-burning phase is similar at Z = 0.020 
and Z = 0.004. Typically, for the 20 M Q models shown on 
Fig. ||, one has that in the non-rotating models, both at 
Z = 0.020 and Z = 0.004, the helium cores contain 26% 
of the total mass at the end of the MS phase. In the ro- 
tating models with V[ n i — 300 km s -1 , this mass fraction 
is enhanced up to values between 30 - 32 %. 

Due to the distribution of the initial velocities and of 
the orientations of the angles of view, rotation induces 
some scatter of the luminosities and effective temperatures 
at the end of the MS phase (see paper V). One observes 
from Fig. H that, at low metallicity, for a given initial ve- 



locity, the extension towards lower effective temperatures 
due to rotation is slightly reduced (compare the tracks for 
«ini = 300 km s _1 ). Thus, at low Z, our models show that, 
if the initial rotation velocities are distributed in the same 
manner as at solar metallicity, the scatter of the effective 
temperatures and luminosities at the end of the MS will 
be reduced. 

6.4. Masses and mass-luminosity relations 

When rotation increases, the actual masses at the end of 
both the MS and the He-burning phases become smaller 
(cf. Tables 1). Typically the quantity of mass lost by stellar 
winds during the MS is enhanced by 45-75% in rotating 
models with v- m i = 300 km s _1 . Similar enhancements were 
found at solar metallicity (see paper V). The increase is 
due mainly to the direct effect of rotation on the mass 
loss rates (cf. Eq. 5). The higher luminosities reached by 
the rotating models and their longer MS lifetimes also 
contribute somewhat to produce smaller final masses. 

Rotation makes the star overluminous for their actual 
masses. Typically for v- m i — 300 km s _1 , the luminos- 
ity vs. mass ratios at the end of the MS are increased 
by 15-22%. It is interesting to mention here that even 
if the rotating and non-rotating tracks in the logL/L Q 
vs. logToff plane are very similar, they may present large 
differences in the log<? c g vs. logT c fj plane where g e s is 
estimated at an average orientation angle as in paper V. 
The large difference in log g c g for similar masses, luminosi- 
ties and effective temperatures comes from the fact that 
the effective gravity of the rotating model differs from the 
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gravity of the non-rotating model by an amount equal 
to the centrifugal acceleration. Typically, the two 40 M Q 
tracks plotted on Fig. |^ which, at \ogT c g = 4.4, differ by 
only ~ 0.04 dex in logL/L©, show important differences in 
the logg c g vs. logT e ff plane. For instance, the rotating 40 
M Q track overlaps the non-rotating 60 M Q model in this 
plane. Therefore, one could expect that the attribution 
of a mass to an observed star position in the log g Q g vs. 
logT ff plane is very rotation dependent. The use of non 
rotating tracks would overestimate the mass (in the exam- 
ple above by 50%), and this might be a cause of the well 
know n pro blem of the mass discrepancy (see e.g. Herrero 
et al. 200C). Let us note that in practice the effective grav- 
ity and the other physical quantities are derived from the 
spectral lines, which shapes and equivalent widths are also 
affected by rotation. 

6.5. Lifetimes 

Table [j] presents some properties of the models. Column 1 
and 2 give the initial mass and the initial velocity v- m i 
respectively. The mean equatorial rotational velocity v 
during the MS phase is indicated in column 3. The H— 
burning lifetimes t-a, the masses M, the equatorial veloc- 
ities v, the helium surface abundance Y s and the surface 
ratios N/C and N/O at the end of the H-burning phase 
and normalized to their initial values are given in columns 
4 to 9. The columns 10 to 13 present some properties of 
the models when logT cr f = 4.0 during the crossing of the 
Hertzsprung-Russel diagram, or when the star enters into 
the WR phase (for the rotating 40 M models and the 
non-rotating 60 M model) , or at the bluest point on the 
blue loop (for the models with M< 12 M Q ). The columns 
14 to 19 present some characteristics of the stellar models 
at the end of the He-burning phase; ijjc is the He-burning 
lifetime. 

From Table one sees that for Z = 0.004 the MS life- 
times are increased by about 7-13% for the mass range 
between 9 and 60 Mq when i>i n i increases from to ^300 
km s _1 . In general, the corresponding changes in the He- 
burning lifetimes are inferior to 10%. As was the case at 
solar metallicity, the ratios tae/ta of the He- to H-burning 
lifetimes are only slightly decreased by rotation and re- 
main around 9-17%. At solar metallicity, the changes of 
the lifetimes due to rotation arc quite similar. The rotating 
20 M Q model with Vi a j = 300 km s^ 1 , at solar metallic- 
ity, has a MS lifetime increased by 14% with respect to 
the non-rotating model. At Z = 0.004, the correspond- 
ing increase amounts to 11%, which is not significantly 
different. 



7. The ratio of blue to red supergiants in the 
SMC 

7.1. Model results 

The observed ratio B/R of blue to red supergiants in the 
SMC cluster NGC 330 lies between 0.5 and 0.8, accord- 



1 i 1 1 1 i 1 — 1 1 i 1 1 1 i 



4.4 - 



4.2 - 



3.8 



3.6 - 



-i 1 1 1 r- 



20 M , Z = 0.004 




0.2 0.4 0.6 0.8 1 

Fraction of the He-burning lifetime 



Fig. 9. Evolution of the T e g as a function of the fraction 
of the lifetime spent in the He-burning phase for 20 M 
stars with different initial velocities. 



i ng to the various sources discussed in Langer & Maeder 
( 1995 ). Not many new results have been obtained since 
then. New IR searches hav e rev ealed some AGB stars 

1996D an d ISO observations 



in the SMC (Zil stra e t al. ^ 

(Kucinskas et al. |2000| ) have led to the detection of an IR 
source in NGC 330, which may be a Be supergiant or a 
post AGB-star, but this does not change the statistics sig- 
nificantly. Notice that the definition of B/ R is not always 
the same, e.g. for Humphreys & McElroy ( 1984 ), B means 
O, B and A-supergiants. Here, we strictly count in the 
B /R ratio the B star models from the end of the MS to 
type B9.5 I, which corresponds to logT c ff = 3.99 accord- 
ing to the calibration by Flower ( 1996 ). We count as red 
supergiants all star models below logT c ff = 3.70 since red 
supergiants in the SMC are not as red as in the Galaxy 
(Humphreys 1979| ). We note that the exact definition of 
this limit has no influence on the observed or theoretical 
B/R ratios, since the evolution through types F,G,K is 
always very fast. 



As noted by Langer & Maeder (1995), the current mod- 
els (without rotation) with Schwarzschild's criterion pre- 
dict no red supergiants in the SMC (cf. Schaller et al. 
1992). This is also seen in Fig. || which illustrates for mod- 
els of 20 M at Z = 0.004 the variations of the T e g as a 
function of the fractional lifetime in the He-burning phase 
for different rotation velocities. For zero rotation, we see 
that the star only moves to the red supergiants at the very 
end of the He-burning phase, so that the B/R ratio, with 
the definitions given above, is B/R ~ 47. For average ro- 
tational velocities v during the MS, v = 152, 229 and 311 
km s _1 , one has respectively B/R = 1.11, 0.43 and 0.28. 
Thus, the B/R ratios are much smaller for higher ini- 
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Table 1. Properties of the stellar models at the end of the H-burning phase, at logT c ff = 4.0 (see text) and at the 
end of the He-burning phase. The masses are in solar mass, the velocities in km s , the lifetimes in million years and 
the helium abundance in mass fraction. The abundance ratios are normalized to their initial value. 
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Fig. 10. Evolution of the T e g as a function of the fraction 
of the lifetime spent in the He-burning phase for 15, 20 
and 25 M stars at Z = 0.004 with u ini = and 300 km 
a" 1 . 

tial rotation velocities, as rotation favours the formation 
of red supergiants and reduce the lifetime in the blue. We 
notice in particular that for v — 200 km s _1 , we have a 
B/R ratio of about 0.6 well corresponding to the range of 
the observed values. 



The B/R ratios change with the stellar masses. Fig. [To| 
shows for the models of 15, 20 and 25 M Q the changes of 
T e ff as a function of the fractional lifetimes in the He- 
burning phase for different rotation. For all masses, we 
notice that the non-rotating stars spend nearly the whole 
of their He-phase as blue supergiants and almost none as 
red supergiants. For i;; n ; = 300 km s -1 (which corresponds 
to about v = 220 km s" 1 ), we notice a drastic decrease of 
the blue phase and a corresponding large increase of the 
red supergiant phase. 

Fig. [ll] shows the same as Fig. ^ but for the models of 
9 and 12 M©. These models mark the transition from the 
behaviour of massive stars, which move at various paces 
from blue to red, to the intermediate mass stars, which 
go directly to the red giant branch and then describe blue 
loops in the HR diagram. At zero rotation, the 15 M Q 
model has the "massive star" behaviour and the 9 M 
model shows a most pronounced "blue loop" . For v ini = 
km s _1 , the 12 M Q model is just in the transition between 
the behaviours of nearby models of 9 and 15 M . The 
rotating model at 15 M is first blue and then goes to 
the red, while the rotating 9 M model goes first to the 
red, then back to the blue and red again. The behaviour 
of the rotating 12 M is also intermediate between these 
two, with the consequence that it always stays more or 
less in the blue, which is surprising at first sight, but well 
consistent with the mentioned intermediate behaviour. As 
seen in Sect. 5, this transition zone with almost entirely 
blue models extends from about 10.5 to 12.2 M . 

The rotating 9 M model has a blue loop smaller than 
for zero rotation; it only extends to the A-type rather than 
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Fig. 11. Evolution of the T e g as a function of the fraction 
of the lifetime spent in the He-burning phase for 9 and 12 
M m stars at Z = 0.004 with u ini = and 300 km s -1 . 



to the B-type range. At a given L and T e ff, the average 
density in a rotating model is much smaller than in the 
non-rotating one, so that the period will be longer. The 
result is such that the application of the standard period- 
luminosity relation will lead for a given observed period to 
a too high luminosity, if the star was a fast rotator on the 
MS. A more complete study of the effects of rotation on 
Ccpheids will be made in a further work. 

Table 2 shows the B/R ratios for the various relevant 
masses for the models with zero-rotation and v- ml — 300 
km s^ 1 . Apart from the transition model of 12 M Q , which 
stays almost entirely in the blue as discussed above, we 
notice that the B/R ratios decrease very much with ro- 
tation, being in the range 0.1 to 0.4 for Vi n i = 300 km 
s _1 . As noted for the 20 M Q model, an average velocity of 
about 200 km s _1 corresponds to a B/R ratio of 0.6. The 
order of magnitude obtained is satisfactory, however, fu- 
ture comparisons in clusters will need detailed convolution 
over the IMF and the distribution of rotational velocities 
in clusters studied at various metallicities. This is beyond 
the scope of this paper and we now examine the effects in 
the internal physics which determine the B/R ratio. 



7.2. Stellar physics and the B/R ratio 

There are several studies on the blue-red motions of stars 
i n the HR diagram, for e xamp le by Lau terbo rn et al. 
(|l97l|) , Stothers and Chin (|l979|) , Maeder ( |l98l|) , Maeder 
fc Me ynet ( 1989 ) and recently by Sugimoto and Fujimoto 
(200C). Sugimoto and Fujimoto identify several param- 
eters at the base of the envelope W, 0, V and N, which 



Table 2. Values of the B/R ratios for models with zero 
rotation and for models with v- ln i — 300 km s _1 . B means 
strictly B-type supergiants and R means K-and M-type 
supergiants (see text). 
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play a role in the redwards evolution. Apart from N, which 
is the polytropic index, we may note that the other pa- 
rameters are all some function of the local gravitational 
potential. V is the ratio of the gravitational potential to 
the thermal energy as in Schwarzschild's textbook (1958). 
The parameter W is given by W = V/U, where U is the 
ratio of the local density to the average internal density. 
The parameter O is given by = ln(— ) c — ln( — ) cnv , where 
the index "c" refers to the center ancf "env" to the base of 
the envelope. We can easily check that is also related to 
the potential at the center and at the base of the envelope, 
as well as to the local polytropic index. 

We may thus wonder whether most of the effects de- 
termining blue vs. red motions in the HR diagram cannot 
be understood, at least qualitatively, in terms of mainly 
the gravitational potential of the core. It is very desirable 
to try to establish some relatively simple scheme for un- 
derstanding the results of numerical computations. The 
role of the core gravitational potential for the inflation 
or deflation of the stellar radius has been emphasized by 
Lauterborn et al. ( 1971 ) in the case of the occurrence of 
blue-lo ops fo r intermediate mass stars (see also Maeder & 
Meynet 198£). We shall examine here whether we may ex- 
tend the very useful "rules" derived by Lauterborn et al. 



(1971) to the case of massive stars in rotation as studied 
here. We call <I> C the potential of the He-core, which due to 
a mass-radius relation for the core behaves as <i>c ~ M® A , 
where M c is the core mass. The blue-red motions in the 
HR diagram mainly depend on the comparison of <!> c with 
some critical potential $ cr it(M), which grows with the 
stellar mass. One has 



$c > *erit(Af) 
$ c < $crit(Af) 



Hayashi line 
blue location 



(8) 
(9) 



Eqs. (8) and (9) essentially apply to a steep hydrogen- 
profile around the core. If this profile is mild, one should 



account for additional terms (Lauterborn et al. 1971). The 
main effect can be represented by a parameter h increasing 
with the amount AA/hc of helium in the transition region 
above the core (formally, a simplified expr ession for h in 
the cases considered by Lauterborn et al. ( |l97l|) is given 
by log/i = 8.5XdMa/M , where is the amplitude of 
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Fig. 12. Comparison of the internal distribution of helium 
in two models of 20 M Q at the middle of the He-burning 
phase. The dashed-dot line concerns the models with zero 
rotation and the continuous line represents the case with 
Wini = 300 km s _1 . 



Fig. 13. Comparison of the internal values of V a( j and of 
V ra d in models of 20 M Q at the very beginning of the He- 
burning phase (Y c — 0.993). The dashed-dot lines concern 
the models with zero rotation and the continuous lines 
represent the case with v- m i = 300 km s _1 . 



the change of the H-mass fraction over the smooth tran- 
sition zone, Md is the width in mass of the transition zone 
and M the total stellar mass). The value of h also depends 
on the distribution of this amount of helium AMh c , h is 
larger if the amount of helium is close to the shell H- 
burning source. However, if a sufficient amount of helium 
is brought far from the shell source, in the outer envelope, 
h may even decrease as suggested by Lauterborn et al. 
(1971). This also has some consequences, as discussed be- 
low. With these remarks, the two above equations become 
for mild H-profiles 



h$ c > $crit(M) 
h<f> c < $crit(M) 



Hayashi line 
blue location 



(10) 
(11) 



We shall now try to see whether we may describe cor- 
rectly with relations (10) and (11) the different physical 
effects influencing the blue-red motions of massive stars: 

Mass loss: Mass loss decreases the total stellar mass 
and thus $ cr i t (M), which favours a motion towards the 
Hayashi line. $ c is not very much changed, since the size 
of the final He-core is not very different. However, there 
is more helium near the H-burning shell, which increases 
the parameter h and also favours the formation of red 
supergiants. 

This description is fully consistent with the well known 
fact that, due to mass loss, the intermediate convective 



zone is much less important (cf. Stothers & Chin |1979 
Maeder 1981 ) . A convective zone imposes a polytropic in- 
dex N ~ 1.5, which implies only a weak density gradient. 



making the stellar radius smaller and thus keeping the star 
in the blue. Thus, the physical connexion we have with 
the interpretation in terms of <i> c is the following one. The 
larger He-burning core with respect to the actual stellar 
mass together with the higher He-content in the H-shcll 
region (higher h) lead to a less efficient H-burning shell, 
thus there is no large intermediate convective zone and 
this absence permits a red location of the star in the HR 
diagram. 

Overshooting: The overshooting does not change 
&ait{M), but <i> c is increased, with no major change of 
the H-profile and thus of h. Clearly, overshooting is thus 
favouring a redwards motion to the Hayashi line, with 
the formation of red supergiants. As for mass loss, the 
larger core contributes to reduce the intermediate convec- 
tive zone, which leads to the formation of red supergiants. 

Lower metallicity Z: A lower Z decreases the mean 
molecular weight (essentially because the He-content is 
also lower, see Sect. 2.3). This decreases the internal tem- 
perature and the luminosity during the MS phase, leading 
to slightly smaller convective cores, as shown by numerical 
models (cf. Meynet et al. 1994). This produces smaller <J> C 
which favours a blue location, as is observed. 

We also note that a lower Z means a slightly higher 
electron scattering opacity (due to the higher H-content), 
which would favour larger cores, however this effect ap- 
pears as a minor one in the models. 

Rotation: The effects of rotation are numerous and 
subtle, and the balance between them is delicate. We no- 
tice the following effects: 



14 



A. Maeder, G. Meynet: Stellar evolution with rotation VII: 




Fig. 14. Comparison of the internal values of V a d and of 
V ra d in models of 20 M Q at the middle of the He-burning 
phase (same model as in Fig. [l2]). Same remarks as in 
previous Figs. 

- 1. A first simple effect is that rotation enhances the 
mass loss rates as described by Eq. (5) and this contributes 
to favour the formation of red supergiants. However, this 
effect although significant does not seem to be the domi- 
nant one. 

- 2. The mixing in the MS phase leads to a slight exten- 
sion of the core, which favours a redwards motion during 
the He-burning phase. This is just like overshooting, thus 
it may be very difficult, except perhaps by asteroseismo- 
logical studies of stars with different rotations, to distin- 
guish between the core extension by overshooting or by 
rotation. An example of this effect can be seen in Fig. |l2|, 
which shows models of a 20 M Q star in the middle of the 
He-burning phase. The core in the non rotating case is 2.7 
Mq, while it is 3.6 M Q in the rotating case, i.e. 1/3 larger 
in mass. Although this increase of the core is very signif- 
icant (in particular for nucleosynthesis), its effect on the 
blue-red motions appears less important than the effect 
we now discuss. 

- 3. A mild mixing just outside the core as produced by 
rotation during the MS phase clearly increases the amount 
of helium near and above the H-shell. This is well seen 
in Fig. |l2|, (in the two models shown, the H-shell is just 
on the right side of the big He-peak). According to the 
definition of h given above, this results in a larger value 
of h, which leads Eq. (10) to be satisfied and this favours 
a red location in the HR diagram. The analysis of the 
sequence of the models before that illustrated in Fig. [l2] 
shows that the different H-profile of the rotating star is 
essentially the consequence of the rotational mixing during 
the MS phase. Thus, the rotational mixing during the MS 



is the key effect for the formation of red supergiants at 
low Z. 

It is satisfactory to see that rules expressed by Eqs. (8) 
to (11) are fully consistent with what we can deduce by 
studying the consequences of a change of mean molecular 
weight. The larger He-core in the rotating models means 
that a larger fraction of the total luminosity is made in 
the core (0.42 instead of 0.31 for the models of Fig. |l2| ). 
This means that the H-burning shell in the rotating model 
produces a smaller fraction of the total luminosity and 
this contributes to reduce the importance of the convec- 
tive zone above the H-burning shell. Simultaneously, the 
higher He content near the shell in the rotating case leads 
to a decrease in the opacity, and this also contributes to 
reduce the importance of the convective zone associated 
with the H-burning shell. Figs. [l3] and |lj show the inter- 
nal values of V ac i and V ra d in 20 Mq stars at the beginning 
of the He-burning phase and at the middle of this phase 
respectively. We notice the smaller convective zone asso- 
ciated with the H-burning shell in the rotating model in 
Fig. |l3] and its earlier disappearance in Fig. |lj. As long 
as it exists, the convective zone maintains a small poly- 
tropic index over the concerned region, and the larger this 
region is, the smaller the radius. Only when the interme- 
diate convective zone disappears, the star reaches the red 
supergiant location in the HR diagram. We see that the 
earlier disappearance of this intermediate convective zone 
in rotating stars favours their earlier evolution towards 
red supergiants, consistent with the result of Eqs. (10) 
and (11). 

- 4. We may wonder what are the respective effects 
of mixing on the MS and during the He-burning phase. 
What we have shown above is just the consequence of the 
mixing during the MS phase. Some tests have been made 
and show the following. If we arbitrarily stop the mix- 
ing during the He-burning phase, this makes little differ- 
ence. However, if on the contrary we arbitrarily enhance 
the mixing during the He-burning phase, this maintains 
the star in a blue location in the HR diagram. This re- 
sult is in agreement with the effect mentioned above and 
discussed by Lauterborn et al. ( 1971 ). The parameter h 
normally grows with AMh c , however it may reach a max- 
imum and then decrease, if helium is brought far enough 
from the H-shcll. Thus, strong mixing in the He-burning 
phase by spreading helium throughout the star may make 
h decrease, which according to Eq. (11) leads to a blue 
location in the HR diagram. 

Why does this occur physically ? Mixing of helium also 
implies mixing of hydrogen. For moderate mixing, the H- 
burning shell becomes more active, due to the increase 
of its H-content and T. This activates the intermediate 
convective zone associated with the H-shell. Additional 
He at the surface also decreases the opacity and favours 
bluer stars. For more extreme cases of mixing both in the 
MS and subsequent phases, one would move to the case of 
almost homogeneous stars, which occurs for fast rotators 
such as the 60 M Q model at Z = 0.02 with a i>i n i of ~ 



400 km s _1 (cf. Meynet & Maeder pOOOj) . In this case, a 
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Fig. 15. Evolution as a function of logT c ff of the abun- 
dance ratio N/C where N and C are the surface abun- 
dances of nitrogen and carbon respectively. The abun- 
dance ratios are normalized to their initial value. The 
tracks are for 20 M Q at Z = 0.004 with different initial 
rotational velocities. The dot-long dashed line shows the 
evolution of a rotating 20 M at solar metallicity (wi„i = 
300 km s" 1 ). 



bluewards evolution in the HR diagram occu rs during the 
MS phase, as shown by Schwarzschild ( |l958|) . 



8. Evolution of the chemical abundances at the 
surface 

As already emphasized in previous works (see Maeder & 
Meyn et 2000a and ref erence s therein; Meynet & Maeder 
200C; Heger & Langer 200C), rotation significantly modi- 
fies the evolution of the surface abundances. Fig. [l| illus- 
trates the changes of the nitrogen to carbon ratios N/C 
from the ZAMS to the red supergiant stage for some 20 
Mq stellar models. For non-rotating stars, the surface en- 
richment in nitrogen only occurs when the star reaches the 
red supergiant phase; there, CNO elements are dredged- 
up by deep convection. 

For rotating stars, N-excesses already occur during the 
MS phase. These N-excesses are not necessarily accompa- 
nied by surface He-enrichments. From Table 1, one sees 
that for the rotating models, well before any change of the 
surface He-abundance occurs, the N/C and N/O ratios 
undergo significant changes. This results from the effects 
of transport on the very strong gradients of the CNO ele- 
ments which rapidly build up inside the star, (we note that 
mass loss at very high rates would do the same). Thus the 
existence of N-rich stars with a normal He abundance is 
easily explained. 



Fig. 16. Surface nitrogen to hydrogen ratios as a func- 
tion of luminosity. The N/H ratios are normalized to their 
initial value. The continuous lines represent evolutionary 
tracks for rotating models with v- m i = 300 km s _1 . The 
initial masses are indicated. The dotted line indicates the 
N/H ratios obtained at the end of the MS phase. The 
dashed line shows this same relation when log T ff is equal 
to 3.9. For the 9 M model which presents a blue loop, 
the value corresponds to the last crossing at logT e ff = 3.9 
before the star becomes a red supergiant. The points rep- 
resent the observed values from Venn (1999), an inverted 
triangle is plotted when only an upper limit is given. 

After the MS phase, the N/C ratios at the surface of 
rotating stars continue to increase during the red super- 
giant phase. At low metallicity, some increase also occurs 
when the star is crossing the HR diagram (see below). The 
N/C ratios obtained at the end of the He-burning phase 
are much higher in rotating models than in non-rotating 
ones. 



From Fig. 15, one sees that at a given metallicity, the 
higher the initial rotational velocity, the more important 
are the surface enrichments at the end of the MS phase. 
This results from the stronger angular velocity gradients 
and transport in faster rotating stars. The increase from 
«ini = 200 to 300 km s _1 produces greater relative changes 
than the increase from 300 to 400 km s~ x (see Fig. |l5| ). 
Likely, this saturation effect results from the compensat- 
ing action of two effects. On the one hand, when the ro- 
tational velocity increases, the transport of the chemical 
species by the shear becomes more efficient. On the other 
hand, when evolution proceeds, this more efficient mixing 
produces smoother chemical and angular momentum gra- 
dients and thus results in a slowing down of the transport 
processes. From Table 1, one sees that, in general, the N 
surface enrichments at the end of the MS phase are higher 
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for higher initial masses. As explained by Maeder (1998) 
this results from the fact that, when the initial mass in- 
creases, the ratio of the mixing timescale for the chemical 
elements to the MS lifetime decreases. 



8.1. Surface enrichments at Z =0.004 and 0.020 

Very interestingly, for a given value of the initial veloc- 
ity, the lower the initial metallicity, the greater the sur- 
face N-enrichments. The mixing of the chemical elements 
is more efficient at low metallicity essentially because the 
metal poor stars are more compact and have greater angu- 
lar velocity gradients. Thus the mixing timescales for the 
chemical elements, which is about Dgbeax/R , are smaller 
at lower metallicities. 

As already briefly mentioned above, the models at low 
metallicity present enhancements of their N/C ratios dur- 
ing their crossing of the HR diagram, while the solar mod- 
els do not present any enrichment during this phase. This 
is a consequence of the fact that at low metallicity, the 
stars begin to burn helium in their core in the blue part of 
the HR diagram. This slows down the evolution towards 
the red and thus gives more time for the transport pro- 
cesses to bring CNO processed material at the surface, 
while the star is still in the blue. 



8.2. Comparisons with the observations 

Fig. |l6| shows the evolution of N/H, the nitrogen to hydro- 
gen ratios, at the surface of our rotating models with 
= 300 km s _1 . The ratios are normalized to their initial 
value. The evolutions of the N/H ratio at the surfaces of 
the non-rotating and rotating stars are qualitatively sim- 
ilar to the behaviours described above for the N/C ratio. 
Let us simply say here that in the plane of Fig. |l6|, the 
non-rotating models predict that stars are either on the 
horizontal line defined by log N/C = or along the line 
log N/C = 0.5-0.6, which corresponds to the position of 
the non-rotating stars having undergone the first dredge- 
up. No stars are predicted outside these two regions. 

Obviously, the observed values obtained by Venn 
( |1999D for A and F supergiants in the SMC are not con- 
centrated along the horizontal lines predicted by the non- 
rotating models. Two observed supergiants show nitrogen 
enhancements well below the values predicted by the first 
dredge-up in non-rotating models. Therefore these can- 
not be on a blue loop. Instead, they are probably on their 
way from the MS to the red giant branch and have un- 
dergone some mixing in the early stage of their evolution. 
A- and F-type supergiants observed by Venn (1999) are 
also observed well above the line corresponding to the first 
dredge-up of non-rotating models. Again, this is an indi- 
cation of an extra mixing process active in massive stars. 

Previous studies have shown that for galactic blue su- 
pergiants, rotation appears as a very promising process to 
account for the observed surface enrichments (see Heger 



& Langer 2000 ; paper V; Meynet 2001 ). Can we say the 



same at low metallicity ? On Fig. [16], the dotted and the 
dashed lines define the evolution between the end of the 
MS phase and the stage corresponding to logT e a = 3.9. 
For stars more massive than about 15 M©, the N-excesses 
do not change very much during this phase. For the stars 
that present a blue loop episode (even a "partial blue loop" 
episode like the 12 M Q model), important enrichments oc- 
cur before the star settles into the red supergiant stage. 

If all the stars were rotating with an initial velocity of 
300 km s _1 , these models would predict that the observed 
A-type supergiants would be around the dashed line in 
Fig. [h| For lower initial velocities the N-enrichments can 
be everywhere between the horizontal line "log N/C = 
" and the dashed line. For higher initial rotational veloci- 
ties, the N-enrichments would be higher than the dashed 
line. Some very N-rich stars at high luminosities (such as 
the N-rich stars with N/H > 10 and log L/Lq > ~ 4.2) 
could be accounted for either by stars with very high initial 
rotation now on their redward tracks, or by very rapid ro- 
tators which, after a red supergiant phase, are going back 
to the blue. Thus we can draw the same conclusions as 
the ones deduced at solar metallicity, namely that there is 
observational evidences for an extra mixing process active 
in massive stars, and that rotational mixing appears as a 
very likely process to drive the extra mixing necessary to 
account for the observations. 

Moreover, the higher enrichments obtained at metal- 
licity lower than solar is in agreement with the observa- 
tions by Venn (1995, 1999). Comparing the range of the 
N/H values measured at the surface of A-type supergiants 
in the Galaxy and in the SMC, she obtained that in the 
SMC, the N/H values cover an interval three times greater 
than in the Galaxy. At logT e ff = 3.9 the N/H ratio at the 
surface of the 20 M© model at Z = 0.004 with v ini = 300 
km s _1 is two times greater than the ratio at the surface 
of the similar model at solar metallicity. For the N/C ra- 
tio the enhancement amounts to a factor 2.5, as can be 
seen from Fig. [l^. Thus, even if the distribution of the ro- 
tational velocities are the same at both metallicities, one 
expects that the low metallicity stars are relatively richer 
in nitrogen. If there are initially more fast rotating stars 
at low metallicity, this will reinforce this trend. 

Thus rotation not only can account for the observa- 
tions of N-rich blue supergiants at both solar and SMC 
metallicity, but it might also account for the fact that 
at low metallicity the maximum enrichments observed are 
greater than at solar metallicity. Finally, we point out that 
in the present models we have found no primary nitrogen 
produced. This is not a difficulty, because the evidence 
for the existence of some primary nitrogen concern metal- 
licities still lower than that of the SMC (cf. Henry et al. 
|2000 ). An interesting question for future models is to see 
whether massive star models with rotation at much lower 
Z produce any primary nitrogen. 

9. Conclusions 

There are 3 main conclusions of this work: 
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-1. The ratios f2/f2 cr i t of the angular velocity to the 
break-up angular velocity grow quite a lot during the 
evolution of massive stars at low Z, contrary to the case 
of massive stars at solar metallicity. This implies that at 
lower Z a larger fraction of massive stars is close to the 
break-up velocity. This effect results from the lower mass 
loss rates and from the significant outwards transport 
of angular momentum by meridional circulation. This is 
very consistent with the observational results by Maeder, 



Grebel and Mermilliod (1999), who observed a larger frac- 
tion of Be stars, when passing from Galaxy, to LMC and 
SMC. The question arises whether the initial distributions 
of the stellar rotation in galaxies of different Z are differ- 
ent. Some arguments in the quoted paper suggest that it 
might be the case, but this is still very uncertain. However, 
we show here that even if the distribution of v sin i is not 
different for massive stars with lower Z, the net result of 
the evolution on the MS leads to much faster rotational 
velocities at the end of the MS phase, when the metallicity 
is lower. 

-2. Due mainly to the additional helium brought near 
the H-burning shell by rotational mixing and the larger 
He-core, which both lead to a less efficient H-burning shell 
and a smaller associated convective zone, the stellar radius 
of rotating stars is permitted to inflate during the He- 
burning phase. These models account for the formation 
of numerous red supergiants at low Z, with a blue to red 
supergiant ratio B/R consistent with the observations in 
the SMC cluster NGC 330. 

-3. While the standard models with mass loss usu- 
ally predict no N/C enrichments on the main sequence 
and in the B- to F-type supergiants, unless they have 
experienced convective dredge-up in red supergiants, the 
models with rotation predict progressive N/C and N/O 
enrichments along the evolutionary tracks. The predicted 
N excesses are in general larger at larger stellar luminosi- 
ties. Also, higher relative excesses of N are predicted at 
lower metallicity by the present models. The comparisons 
with the observations by Venn (1998, 1999] ) of A-type su- 
pergiants show a good agreement for the size of the N- 
enrichments, as well as for the fact that they are larger at 
lower Z. 

The major interrogation for the future lies in the line of 
the first conclusion. Is rotation more and more important 
as we go to massive stars of lower Z ? This is a question 
for which great observational efforts need to be made. If 
the answer is yes, for whatever reason (star formation or 
effect of MS evolution), this will imply that rotation is a 
dominant effect and has important consequences for the 
interpretation of the spectra of high redshift galaxies and 
for the early nucleosynthesis. 
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